Electron shot noise beyond the second moment 
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The form of electron counting statistics of the tunneling current noise in a generic many-body 
interacting electron system is obtained. The third correlator of current fluctuations (the skewness 
of the charge counting distribution) has a universal relation with the current I and the quasiparticle 
charge e*. This relation C3 = (e*)^/ holds in a wide bias voltage range, both at large and small 
eV/kB,T, thereby representing an advantage compared to the Schottky formula. We consider the 
possibility of using the counting statistics for detecting quasiparticle charge at high temperature. 
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Recent developments in the problem of quantum elec- 
tron transport were marked by interest in the phe- 
nomenon of electric noise. The many-body theory of elec- 
tron shot noise, developed by Lesovik Q| (and indepen- 
dently by Khlus [^) for a point contact, was extended 
to multiterminal systems by Biittiker |^ and to meso- 
scopic systems by Beenakker and Biittiker Q . Kane and 
Fisher proposed using shot noise for detecting fractional 
quasiparticles in a Quantum Hall Luttinger liquid |^]. 

Experimental studies of the shot noise, after first mea- 
surements in a point contact by Reznikov et al. j^] and 
Kumar et al. [Q, focused on the quantum Hall regime. 
The fractional charges e/3 and e/5 were observed H-[lo|] 
at incompressible Landau level filling (see also recent 
work on noise at intermediate filling ) . The shot noise 
in a mesoscopic conductor was observed by Steinbach et 
al. IT2) and Schoelkopf et al. |T^, who also studied noise 
in an ac driven phase-coherent mesoscopic conductor . 

In this article we discuss a generalization of the shot 
noise, namely the counting statistics of fluctuating elec- 
tric current. It can be defined through the probability 
distribution P{q) of charge transmitted in a fixed time 
interval ||l^,|f^ . We consider ways of obtaining the distri- 
bution P{q) using a fast charge integrator scheme. From 
the distribution P{q) all moments of charge fluctuations 
can be calculated and, conversely, the knowledge of all 
moments is in principle sufficient for reconstruction of the 
full distribution. However, due to the central limit theo- 
rem, high moments are difficult to access experimentally. 
Therefore we shall focus primarily on the third moment. 

The counting statistics have been analyzed theoreti- 
cally for a Fermi gas, in the single- and multi-channel 
geometry ]l5| , p7[ , in the mesoscopic regime [p^p^ , and 
in the ac driven phase-coherent regime ||r^ , |20|| . Charge 
doubling due to Andreev scattering in NS junctions was 
considered by Muzykantskii and Khmelnitskii |2^, and 
in mesoscopic NS systems by Belzig and Nazarov [ p2[ . 
However, since the most interesting applications of the 
shot noise lie in the domain of interacting electron sys- 
tems, an appropriate extension of the theory is necessary. 

The problem of back influence of a charge detector 
on current fluctuations was considered by Lesovik and 



Loosen p^ , and recently by Nazarov and Kindermann 
[p^ . Beenakker proposed an alternative way of obtaining 
charge statistics using photon counting [ p5[ . Application 
to pumping in quantum dots was also discussed . 

Our central finding is a relation between counting 
statistics and the Kubo theorem, valid in the tunneling 
regime for a generic interacting many-body system. From 
that we obtain a formula for the moments of the counting 
statistics that holds in the entire bias voltage range, at 
arbitrary eV/k^T. We demonstrate that in the tunnel- 
ing regime the current fluctuations are described by an 
uncorrelated mixture of two Poisson processes. This is re- 
vealed by a generating function x(A) = P(q)e''*"^/'^ , 
with e* the quasiparticle charge. We find 



X(A) =exp [(e»^-f)A^i^2(r) + (e 



-l)N2 



0] (1) 



where Na^bij) = niabT is the mean charge number trans- 
mitted from the contact a to the contact 6 in a time r. 

The result (|I|) yields a number of relations between 
different statistics of the probability distribution P{q). 
The cummulants {{Sq'^)) (irreducible correlators) of the 
distribution P{q) are expressed in terms of xW as 



Inx(A) = 



(»A)^ {{Sq'^)) 
k\ {e*Y 



Using Eq. (hi) one obtains 



' [ [mi2 + m2i)r. 



(2) 



k odd 
k even 



Setting A; = 1, 2 we express TO12 ± m2i through the time- 
averaged current and the low frequency noise power: 



*\2 



TO12 - m2i = //e*, mi2 -I- m2i = 7^/2(6*) 



(4) 



Of special interest for us will be the cummulant {{Sq'^)) 
which is equal to the third correlator |2^] 



_n3 



{{6q'))^6q^ = {q-q) 



(5) 
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(see Fig. ||). For this correlator Eq. gives {{Sq^)) = 
C3T with the coefficient C3 ("spectral power") related to 
the current / as 



(6) 



We note that the relation ^ holds for the distribu- 
tion (|^) at any ratio of the mean number of transmitted 
charges mi2 — 11121 to the variance mi2 + m2i. 
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FIG. 1. The third moment (^) determines the shape of the 
distribution P{q), namely its skewness. This is illustrated by 
a distribution of the form (0) and a Gaussian with the same 
mean and variance. For C3 > the peak is somewhat more 
stretched to the right than to the left. 

The meaning of Eq. ^ is similar to that of the Schot- 
tky formula for the second correlator V = 2{{5q^J) = 2e*/ 
which is usually used to determine the effective charge e* 
from the tunneling current noise. The Schottky formula 
is valid when charge flow is unidirectional, which means 
mi2 ^ TO21 (see Eq. (||)). The latter can be true only at 
sufficiently low temperatures fcsT' <^ eV. This require- 
ment of a cold sample at a relatively high bias voltage is 
the origin of a well known difficulty in the noise measure- 
ment. In contrast, the relation (^) is not constrained by 
any requirement on sample temperature. 

On a general basis we expect the relation (^ to hold 
approximately even outside the tunneling regime. In- 
deed, for the Nyquist noise at equilibrium all odd mo- 
ments vanish. Combined with the temperature indepen- 
dence of P{q) out of equilibrium, this implies a relatively 
weaker dependence on eV/kBT than in the noise power 
at the Nyquist-Schottky crossover. This is manifest, for 
instance, in the temperature independent first moment 
of P{q) for free fermions (the Landauer formula). 

This property of the third moment, if confirmed exper- 
imentally, may prove to be quite useful for determining 
the quasiparticle charge. In particular, this applies to 
the situations when the I — V characteristic is strongly 
nonlinear, when it is usually difficult to unambiguously 
interpret the variation of the second moment with cur- 
rent as a shot noise effect or as a result of thermal noise 
modified by non-linear conductance. We stress that this 



is a completely general problem pertinent to any interact- 
ing system. Namely, in systems such as Luttinger liquids, 
the I — V nonlinearities arise at eV > kBT. However, it is 
exactly this voltage that has to be applied for measuring 
the shot noise in the Schottky regime. 

Now we turn to the derivation of the main result (|^). 
The starting point of our analysis will be the tunneling 
Hamiltonian H — Hi + 'H2 + V, where 7ii^2 describe 
the leads and V — J12 + J21 is the tunneling operator. 
The specific form of the operators J12, J21 that describe 
tunneling of a quasiparticle between the leads will be 
unimportant for the most of our discussion. 

The counting statistics generating function xW can 
be written pq| as a Keldysh partition function 



X(A) = (Tk 



exp 



nx{t)dt 



(7) 
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where a counting field X{t) is added to the phase of the 
tunneling operators J12, J21 as 



Vx = ei^(*)Ji2(t) + e-5A(t)^(i) 



(8) 



Here X(t) — ±A is antisymmetric on the forward 
and backward parts of the Keldysh contour Co,t = 
[0 ^ r ^ 0]. Eqs. (|^), (||) originate from the analysis 
of a coupling Hamiltonian for an ideal "passive charge 
detector" without internal dynamics p6|j2^ ]. 

In what follows we compute x(A) and establish a re- 
lation with the Kubo theorem for tunneling current [ p9[ . 
For that, we perform the usual gauge transformation 
turning the bias voltage into the tunneling operator phase 
factor as J12 — > ^126"'^^^*, J21 J2ie*'^^'. Passing to 
the Keldysh interaction representation, we write 



(9) 



X(A) = (Tkcxp -I / Vm)it)dt 



Diagrammatically, the partition function (^ is a sum of 
linked cluster diagrams with appropriate combinatorial 
factors. To the lowest order in the tunneling operators 
J\i, J21 we only need to consider linked clusters of the 
second order. This gives x(A) = e^^'^-', where 



W{X) = -- 



TKVxit){t)Vxit'){t')) dtdt' (10) 



There are several different contributions to this integral, 
from t and t' on the forward or backward parts of the 
contour Co.r- Evaluating them separately, we obtain 



W{X) 



V-x{t)Vx{t'))dt'dt 



(11) 



yx{t)Vxit')) dt'dt - J (V-xit)V-x{t')) dt'dt 



We substitute the form (^ into Eq.(ll) and average by 
pairing J12 with J21. This gives 
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with N„^h 



2(T) + (e-^^-l)iV2_i(r) 

r 

{Jba{t)Jab{t'))dtdt' 



(12) 
(13) 



Exponentiating (|T^) gives the result ([^) 

It is instructive to relate the quantities (13) with the 
Kubo theorem. We consider the tunneling current oper- 
ator I{t) = -ie* (^Ji2it) — J2i(i)^- From the Kubo the- 
orem for the tunneling current |p9[ , the mean integrated 
current {X{t))dt scaled by e* is nothing but 



■hi{t)Ji2{t') ) dtdt' ^ Ni^2 - N2^l (14) 



^0 

By writing Na^b ~ fn-abT, we obtain the first relation 
(Q). To obtain the second relation (^ we consider the 
variance of the charge transmitted in time r. It is given 
by a time integral of an averaged symmetrized product 
of two current operators PSI 



{{Sq')) = (e 



*\2 



|i(i),Z(t')} )dtdt' (15) 



The integral in (|T^) can be rewritten as 

^ ^ {{Jl2{t), J2l{t')}^) dtdt' = Ni^2+N2 



(16) 



which immediately leads to (^). 

The result (^, and thereby the formula (^ for the 3-rd 
correlator, are valid only at low transmission. In that the 
situation is similar to the Kubo formula for the tunneling 
current, which is valid only in the tunneling Hamiltonian 
approximation. To illustrate this we recall the expression 
for counting statistics for a single channel noninteracting 
Fermi system (point contact]in the presence of a dc volt- 
age V and temperature T [fa, 



X(A) = exp(-iVTW+W_) 



2'Kh 



(17) 



U± = U/2 ± cosh"^ (t cosh(W/2 +r cosh {U/2)) (18) 

where r is a measurement time, and U — eV/k^T . This 
result holds for any values of the transmission and reflec- 
tion constants t and r (constrained by i -I- r = 1). The 
formula (|l^) was obtained in Ref. |l^ by explicitly eval- 
uating the Keldysh partition function in the scattering 
basis representation. 

The 3-rd correlator {{Sq^)) can be obtained from ( |l^ ) 
by expanding Inx(A) in Taylor series up to O(A^): 



{{Sq-")) ^ e%l - t)NT [ 6t 



sinhW — U 

L 

coshZ// — 1 



+ {l-2t)U) (19) 



This expression is a function of the bias-to-temperature 
ratio hi, and so in this case the relation ( p^ ) for the 3-rd 
correlator does not hold (see Fig. ||). Asymptotically 



where / 



e2(l 
e2(l 



t)lT, 

20(1 



t)/r. 



eV < kBT 
eV > kBT 



(20) 



T^tV. One can also average over the uni- 
versal Dorokhov's distribution of transmission in a mul- 
tichannel mesoscopic metal |l|] ( see Fig. Eh. 
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FIG. 2. The third correlator C3 scaled by e /, with / the 
time-averaged current [see Eqs. (^), (^], for the single chan- 
nel problem (p^), (p^, and for a phase-coherent mesoscopic 
multichannel conductor. Note that the relation C3 = e^I 
holds approximately at not too large transmission t. 

Eqs. (p^, ( |2^ ) indicate nonuniversality of the relation 
(^) outside the tunneling regime. They also lead to an 
interesting qualitative prediction: At t > 0.5 the ratio 
{{5q^J)/I can become negative. Such a signature could 
be observed even if it proves difficult to measure {{Sq^)) 
quantitatively with sufficient precision. This is impor- 
tant in view of the difficulties in measuring the counting 
statistics (see below). 

In the single channel problem ( p7| ) the tunneling regime 
is realized at low transmission t. To connect with the re- 
sults (|I]), (||) we analyze the expression (|l^) at t -C 1. To 
the lowest order in small t we have 



e'^ - 1 



-1 



U+ ^U, U-^ t- 

- i 

Substituting this in Eq. (p^) we recover (Q) with 
eVr t 



(21) 



N2^liT) 



2Trh -1' 



iVi^2(T) =e"iV2^i(T) (22) 



the rates of two Poisson processes. 

The measurement of the distribution P{q) is a nontriv- 
ial task. Current fluctuations must be amplified with a 
very low noise preamplifier {e.g. the one used in Refs. 
[^,|o|). Amplified signal can then be digitized and ana- 
lyzed with computer. This setup in principle allows to 
reconstruct the full statistics of transmitted charge. In 
practice, however, the correlators of high order become 
increasingly difficult to extract. 

The main source of error in the measurement of the k- 
th cummulant Ck of the distribution P{q) is statistical. 
The nongaussian character of the amplifier noise does not 



3 



present a problem, since the mean time-averaged value 
of the fc-th cummulant c|,'^^ for the amplifier can be sub- 
tracted if known with sufficient accuracy. The measured 
value Cfc should be compared to (i) the variance var Ck 
of the fc-th cummulant statistics and (ii) the variance 
var C^"'' of the fc-th cummulant of the amplifier noise. 
The variance is in both cases expressed through the cor- 
relators of order 2k. The correlators of even order for a 
generic distribution can be estimated, by virtue of the 
central limit theorem, using Gaussian statistics: 



var Ck - {{Sq^"'))'^^ ^ ((2fc - 1)!! {Sq')") 



1/2 



(23) 



Fluctuations introduced by amplifier can also be esti- 
mated using Gaussian statistics. For the amplifier noise 
of spectral density A (measured in A^/Hz), charge fluctu- 
ations are SQ^ = ^At, where t is sampling time. An es- 
timate of the variance var C^"', similar to Eq.(|23|), gives 



var C, 



(a) 



((2fc-l)! 



'Mr) 



1/2 



(24) 



For odd k the commulant Ck — {{^q'^)) mean value is 
proportional to the current /, as discussed above. The 
fluctuations due to the amplifier are independent of /. 
In the Nyquist regime (at small /) the variance var Ck is 
independent on / and is determined by thermal noise. In 
the Schottky noise regime (at large /) the fluctuations are 
determined by /, and var Ck oc J*^/^. Therefore at fc > 2 
one can increase the signal-to-noise ratio by increasing 
the current until var Ck ~ var C^°^ 

In the shot noise regime, when {dq'^) = elr 
can estimate the signal-to-noise ratio as 



At, we 



S/N: 



1 



2varCfc 2((2fc - l)!!)i/2 



2\ fc/2-1 



2e2\ 
A^ 



(25) 



It is clear from Eq. ( |25| ) that it is beneficial to decrease the 
sampling time r to gain sensitivity. Estimates for typical 
values of A = 10"^^ A^/Hz and r = 10^'' s~^ give for 
the fc = 3 commulant S/N « 10^^. This value is accept- 
able, since repeating the measurement many times over a 
long time T and averaging w ill fur ther reduce statistical 
fluctuations by a factor of -y/T/r. For the cummulants 
of higher order fc > 3 the situation is more problematic. 

In summary, the counting statistics (|l|) of tunneling 
current is found to be universal and independent of the 
character of interactions. For the third correlator we ob- 
tain a generalized Schottky formula (^. This formula is 
valid at both large and small eV/kBT and can be used to 
measure quasiparticle charge at temperatures fceT > eV. 
A method for measuring the third correlator is proposed. 
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